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ABSTRACT 

The Magnetodynamics, or Force-Free Degenerate Electrodynamics, is recognised as 
a very useful approximation in studies of magnetospheres of relativistic stars. In this 
paper we discuss various forms of the Magnetodynamic equations which can be used 
to study magnetospheres of black holes. In particular, we focus on the 3+1 equations 
which allow for curved and dynamic spacetime. 
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1 INTRODUCTION 

In magnetospheres of pulsars and black holes the electro- 
magnetic field is so strong that inertia and pressure of 
plasma can be ignored. As the result, the Lorentz force 
almost vanishes and the transport of energy and momen- 
tum i s almost entirely elec t roma gnetic (|Goldreich fc Julian! 
1 19691 : iBlandford fc Znaiekl [l977l ). This justifies the name 
"force-free" to describe the electrodynamics of pulsars and 
black holes. However, the electrodynamics of the magneto- 
spheres is rather different from the electrodynamics in vac- 
uum which, obviously, is also force-free. Indeed, the mag- 
netospheric plasma is plentiful enough to support strong 
electric currents and screen the component of electric field 
parallel to the magnetic field. Electromagnetic field satisfy- 
ing this condition is ca l led "d egenerate" and for this reason 
iMacdonald fc Thorne I called the electrodynamics of 

pulsar and black hole magnetospheres "force-free degenerate 
electrodynamics" (FFDE). 

For a long time theorists were preoccupied with steady- 
state solutions of FFDE. Even the basic properties of FFDE 
as a system of time-dependent equation were not studied 
systematically. The first ste p in t his di rection wan made 
only quite recently, when lUch idal |l993) developed a the- 
ory of FFDE in which the electromagnetic field is described 
in terms of two scalar functions, called "Euler potentials". 
However, this formulation has not been very popular. In 
particular, it is not very convenient for numerical analysis 
because its basic equations, when written in components, 
involve mixed space and time second order derivatives. An- 
other approach is to use the actual Ma:xwell equations sup- 
plemented with a parti cular pres c riptio n for the electric cur- 
rent. This was done bv lGruzinovl (|l999l ), who used th e force - 
free condition to derive the Ohm law. iKomissarovl (|2002l ') 
showed that FFDE can be considered as Relativistic MHD 
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(RMHD) in the limit of vanishing inertia of plasma parti- 
cles. This allowed to rewrite FFDE as a system of conserva- 
tion laws similar t o RMHD, including the en ergy-momentum 
conservation law. iKomissarov et all (|2007l l argued that the 
dynamics of electromagnetic field in FFDE can be inter- 
preted as a motion of magnetic mass-energy under the action 
of Maxwell stresses and proposed another name, "Magneto- 
dynamics" (MD), for FFDE. We will be using this name in 
the rest of the paper. 

The formulation by IKomissarovl (|2002l ') is in a co- 
variant form and can be used to s t udy t he magne- 
tosph eres of black holes (|Komissarovl l200ll : iMcKinnevI 
l2006f ). However, the wealth of experience accumulated 
in solving Maxwell equation s has ensured that the for- 
mulation by GruzinovT ( 19991) was often found preferabl e 
l|Spitkovskvl l2006l : iKalapotharakos fc ContopoulosI 120091 '). 
This prompted recent efforts to generalise Gruzinov's for- 
mulation so that it could also be used to study the magne- 
tos pheres of black holes. The s tarting point was the work 
bv lThorne fc Macdonaldl (| 19821 ') . who first obtained general 
3+1 equations of Electrodynamics (Eqs.(3.4) of this pa- 
per) and then a simplified version (Eqs.(5.8) of the same 
paper) which was adapted to the case of stationary black 
holes. Moreover, they restricted their attention to the Boyer- 
Lindquist foliation of spacetime. This simplified version has 
become most known to astroph y sicists , via the follow-up 
paper by IMacdonald fc Thornel (Il982l) a nd iThorne et"al] 
(| 19861 ). and widely used. Komissarovl (|20o4 ) develope d a dif- 
ferent formulation, which has its roots in the works o f'Tam: 
( 1924) and Pleb anski ( 1959 ) (see also L andau fc Lif shit 
(119711 )). In this formulations the 3+1 equations of Electro- 
dynamics also have a very simple and familiar form. In fact, 
they look exactly the same as the Maxwell equations in mat- 
ter. The only assumption on the spacetime metric made in 
this formulation is that the determi nant of the metric tensor 
of space does not depend on time. iPalenzuela et all (|201G| ) 
presented, without derivation, the 3+1 equations which are 
free even from this constraint. They seem to have used the 
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approach bv lKomissarov (|2004) but rev e rted t o the original 
representation of Thorne fc Macdonaldl l|l982l ). where only 
the electric E and magnetic field B are present. Their equa- 
tions also include extra scalar fields, which have been intro- 
duced for purely computational reasons. The force-free Ohm 
law of General Relativistic Magnetodyna mics (GRMD) was 
first derived in the space-time f orm by McKinneyl (|2006| ) 
and then in t he 3-1-1 form by IPalenziiela et alj (|201lh . 
iLvutikovl (|201ll ) independently derived equations of GRMD 
using the simplified version of 3-1-1 Electrodynamics by 
[rii orne &: Macdonaldl (|l982l ). Thus, his equations have in- 
herite d the limitations of those by iThorne fc Macdonaldl 
j 19821 ). 

In this paper, we revert back to the 3-1-1 formulation 
of I Komissarov! (|2004l '). modify it in order to allow non- 
stationary metric, and derive the corresponding form of 
the force-free Ohm law. We also present various relevant 
derivations and explore the connections between the differ- 
ent forms of the 3-1-1 equations. 



2 3+1 ELECTRODYNAMICS 

Following iThorne fc Macdo nald' (1982') we adopt the folia- 
tion approach to the 3-1-1 splitting of spacetime in which the 
time coordinate t parametrises a suitable filling of spacetime 
with space-like hypersurfaces described by the 3-dimensional 
metric tensor 7ij. These hypersurfaces may be regarded as 
the "absolute space" at different instances of time t. Below 
we describe a number of useful results for further references. 
If {x^} are the spatial coordinates of the absolute space then 



ds 



(/3^ — a^)dt^ + 2Pidx^dt + fijdx'dx-' , 



(1) 



where a is called the "lapse function" and (3 is the "shift 
vector" . The 4-velocity of the local fiducial observer, 'FIDO', 
which can be described as being at rest in the absolute space, 
is 



= (-Q, 0,0,0). 



(2) 



The spatial components of the projection tensor, which is 
used to construct pure spatial tensors, 



(3) 



coincide with the components of the spatial metric . Other 
useful results are 



a 



g = 

a 



g ^ -a J , 
where 



(4) 

(5) 
(6) 



P' ^f^Pj, g = detg^^, 7 = det 7,^ 



are the components of the velocity of the spatial grid rela- 
tive to the local FIDO as mea sured using the coordinate tim e 
t and the spatial basis {ft} (|Macdonald fc Thorne I [l98l ). 

The covariant Maxwell equations are (e.g. Jack- 
son(1975)): 

V;3*F'^''=0, (7) 
and 



VpF"^ = r , (8) 

where F"'^ is the Maxwell tensor of the electromagnetic field, 
*pai3 ^jjg Faraday tensor and I"' is the 4-vector of the 
electric current. The most direct way of 3+1 splitting of 
the covariant Maxwell equations is to write them down in 
components and then to introduce such spatial vectors that 
these equations have a particularly simple and familiar form. 
For example, when Eq.(l7]) is written in components it splits 
into two parts: 



The time part: 

d,(^a^*F") =0, 



^/7 

• The spatial part; 
-^ft + -^ft (aV7*f^') =0. 

If we now introduce the spatial vectors B and E via 



-e-^F^'' 



and 
E, = 
where 



ijfe 



(9) 

(10) 

(11) 
(12) 

(13) 



is the Levi-Civita tensor of the absolute space and e^jfe = 
e'-'*' is the 3 dimensional Levi-Civita symbol, then equations 
(f9lT0ll read 

V-B = 0, 

1 



^/7 



9t(y7B) + VxB^O, 



(14) 
(15) 



where V is the covariant derivative of the absolute space. 
Similarly, equation ^ splits into 



^9t(y7£')- VxJf = -J, 
where 

and 

p — OiL ^ J — ai . 



(16) 
(17) 

(18) 
(19) 

(20) 



Similar to any highly ionised plasma, the pair plasma of 
black hole magnetospheres has essentially zero electric and 
magnetic susceptibilities. In such a case, the Faraday tensor 
is simply dual to the Maxwell tensor 



F^v 



Here 



(21) 
(22) 

(23) 
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is the Levi-Civita alternating tensor of spacetime and 
is the four-dimensional Levi-Civita symbol. 



This allows us to obtain the following alternative expressions 
for B,E,D, and H: 



D' = ie 



ijk 



jk , 



(24) 
(25) 
(26) 
(27) 



Moreover, from the above definitions one immediately finds 
the following vacuum constitutive equations: 



E ^ aD + l3xB , 
H = aB pxD. 



(28) 
(29) 



In flat spacetime with Lorentzian (pseudo-Cartesian) coor- 
dinates one has a — 1, P = and, hence, B — H and 
E = D. 

Each of the introduced spacial vectors can be repre- 
sented by a spacetime vector whose spacial part is the spa- 
cial vector in question and whose time part vanishes. As one 
can easily verify, these spacetime vectors are given by the 
following covariant expressions: 



f 



(30) 

(31) 
(32) 

(33) 

(34) 



where fc" = dt- From these one can see that B and D are the 
magnetic and electric fields as measured by FIDOs, whereas 
H and E are auxiliary vector fields. 
It is also easy to verify that 



(35) 



and thus p is the electric charge density as measured by FI- 
DOs. However, J, is not the electric current density as mea- 
sured by FIDO, which we will denote as j. Geometrically, j 
is the component of I" normal to n". Using the projection 
tensor 71, ^ = g,y^ + n,,n^, we find 



J ^ aj - pf3 . 



(36) 



The second term in this equation accounts for the motion 
of spatial grid relative to FIDO, or in other words for the 
fact that the coordinate time direction, the basis vector dt, 
is generally not parallel to n". 

When dtj = these 3-1-1 equations have exactly the 
same form as the classical Maxwell equations for the elec- 
tromagnetic field in matter 



V-B = 0, 

dtB + VxE = 0, 

V D ^ p, 

- dtD + VxH ^ J . 



(37) 
(38) 
(39) 
(40) 



This similarity explains why we prefer to denote the electric 
field measured by FIDO as D, whereas in most papers by 
other researches it is denoted as E. 

Applying the divergence operator to Ea lI7l one finds the 
electric charge conservation law 



I 



V7 



9t(%/7p) + V-J = 0. 



(41) 



Although this is slightly different from the usual differential 
form of this law, its integral form is exactly the same 



J dS = 0, 



(42) 



where dV is the metric volume and dS is the metric surface 
elements. 

The limit of Magnetodynamics is defined by vanishing 
of the Lorentz force. In the covariant form this condition 
reads as 

Fm./^ = . (43) 
In our 3+1 formulation this equation splits into 
E-J^O (44) 
and 

pE + JxB = 0. (45) 
From the last equation it follows that 

E-B^O. (46) 



When combined with the constitutive equation (|28|l , the last 
equation also implies 



D B = 0. 



(47) 



As first noticed bv lCruzinovl (|l999l ). the force-free condition 
allows one to express the electric current in terms of the 
electromagnetic field and its spacial derivatives, thus pro- 
viding us with a particular form of Ohm's law. Here we re- 
peat Gruzinov's derivation taking into account the effects of 
General Relativity. The component of electric current nor- 
mal to the magnetic field can be found directly from Eg 1451 
via cross-multiplying its sides by B. This yields 

= P^ . (48) 

In order to find the parallel component we first notice that 
EaliTl implies 

dt{VlD-B) = 0. (49) 

When combined with Eqs. (|I5ll7l47|) this yields 

{VxH - J)- B - (VxE) ■ D = 0, (50) 

which does not involve the time derivative of 7. From the 
last result we find that 

j^^^ B.iVxH)^D.i^xE) ^^ (51) 

Collecting all these results, we can write the most gen- 
eral 3-1-1 system of GRMD as 



V7 



dti^B) + VxE^O, 



(52) 
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—dt{y^D)-VxH = -J, 

V B = 0, 

where 

E = aD + f3xB , 

H = aB ~ /3xD , 

ExB B ■ CVxH) - D ■ {VxE) ^ 
J ~ P — — 1 B , 



£2 



£2 



and 



(53) 
(54) 

(55) 
(56) 

(57) 
(58) 



r) = aTr(K), 



(63) 



It is easy to see that in flat spacetime with Lorentzian 
coordinates, where q = 1, /3 = 0, an d dt^ = 0, this system 
is reduced to that of lGruzinovl l|l999l ). Under the c onditions 
9f7 = and V • /3 = it is reduced to that of iLvutikovl 
ll201iri . 

Magnetodynamics can be considered as Relativistic 
Magn etohydrodynamics in the limit of vanishing particle in- 
ertia CKomissarov 2002). The explicit condition of magne- 
tohydrodynamic approximation is vanishing of the electric 
field in the fluid frame. This implies that in any other frame 
the component of electric field parallel to the magnetic one 
always vanishes and the magnetic field is stronger than the 
electric one. These conditions can be written in the covariant 
form as 



*F^^F'"' = and 



> 0. 



In our 3+1 notation these yield 
B D = and - > . 



(59) 



(60) 



In computer simulation, one has to make sure that the ini- 
tial solution satisfies both these conditions. The first con- 
straint is preserved exactly by the differential equations 
of MP. However, the second constraint can be violated 
l|Komissarovl | 2002! ). Slow shocks of RMHD can transform 
plasm a from magnetic ally-dominated to particle-dominated 
state l|Lvubarskvi2005l ). Ho wever, slow waves a re not allowed 
in the MD approximation (lKomissarovll2002l ). This limita- 
tion can be behind many violations of the second condition 
in MD. 

Substituting the expressions for E and H from the con- 
stitutive equations into Eas. (|15ll7|l and expanding the dou- 
ble cross-product terms one finds 



dtB - CpB + Vx aD = rjB 

and 

dtD - jC/sD -VxaB^rjD-aj 



(61) 



(62) 



where rj = ^ ■ f3 — 9t(ln^/7) and C/s is the Lie derivative 
along the shift vector (e.g. CfsB = (/3-V)B - (B- V)/3 
). This is another useful form of the most general Fara- 
day and Ampere equa tions of 3+1 GR Electrodynamics 
l|Palenzuela et ahlboid fl. One can show that 



^ In | PaIenzueIa et al. as well as in lThorne fc Macdonaldl 

and many other papers, the variable D is denoted as E, 
following its interpretation as the electric field measured by the 
local FIDO of the spacetime foliation. 



where Tr{K) = 'y'''°Kik is the trace of the external curvature 
tensor of the absolute space 



Kik = — {Pi-k + Pk;i — dtjik) 

2a 



(64) 



(|Misner et al.ll 19731 '). When both dt^ = and V-/3 = these 
equations reduce to 



dtB - CpB + V X aD = 

and 

- dtD + CfiD + VxaB = aj . 



(65) 



(66) 



These are the 3+1 equations of Bla ck Hole Electrodynam- 
ics bv iMacdonald fc Thorne I (|l982l '). We note here that al- 
though the condition V • /3 = is satisfied by the Boyer- 
Lindquist metric of Kerr black holes it is not satisfied by 
the Kerr-Schild metric, which is also widely used in black 
hole studies. 

In terms of the physical quantities measured by FIDOs, 
vanishing of the Lorentz force has the familiar fornfl 



pD + jxB ^ . 



(67) 



The force-free electric current j can now be obtained in ex- 
actly the same fashion we did earlier for J. The normal 
component of j is obviously 

DxB , , 

= P^i- ■ (68) 

In order to find the parallel component we apply the oper- 
ator dt — Cp to B-D = 0. This yields 



3ll = 



B ■ ( V X aB) - D ■ (qV x D) 



(69) 



aB2 

Given the identity A ■ (V x aA) = aA ■ (V x A) the final 
expression for the force free current does not actually involve 
either the shift vector or the lapse function and has exactly 
the same form as in Special Relativity, 



. _ DxB B- (VxB) - £>■ (VxD) 

3 - P ^2 + 52 

iPalenzuela et al.ll20nT ). 



(70) 



3 THE 4- VECTOR OF FORCE-FREE 
CURRENT 

Finally, we briefly discuss the space-time formulation of 
Magnetodynamics. If one prefers to deal with the 4- 
tensor Maxwell- Ampere e quation jSl) instead of the energy- 
momentum equation (as in lKomissarovl (|2002l ll then the key 
issue is the expression for the 4- vector of force-free current. 
This expression was found by iMcKinnevI (jiboQ l. However, 
it can be simplified a little bit further. Here we explain this 
and give a slightly different derivation. 

From the definitions (I3QI32[) it follows that 



(71) 



and 



2 This equation can also be obtained via substituting expressions 
ll28l l and (|36]l into Eg.BSt. 
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-n B' + B n ~ e Di,n^ . 



Then the force-free condition (|43|l reads 
pD^ + e^^^^'nJ^B, = 0. 



From this we find that 
3/3 



Z^^Sfl = 



and 

B^ 



n-,DpBs + ^^-^B^ + pn" 



(72) 



(73) 



(74) 



(75) 



One can see that the spacial part of which we will denote 
as J'^ = J^+J^i has the following components parallel and 
perpendicular to 



B2 



(76) 
(77) 



The coefficient I'^B^ in Ea l76l can be expressed in terms 
of the electric and magnetic fields and their derivatives, mak- 
ing this equation an explicit expression for Following 
iMcKinnevI i 20061 ) we first contract the Maxwell- Ampere law 
(HI with to find that 



rB^ = -B"D. 



(78) 



where the comma indicates partial derivative. Then we con- 
tract the Maxwell-Faraday equation ([7]) with to find that 



Thus, 



OLdfB 7-1 7-1 



(79) 
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7" Be = e^°'^-'n^{Bc,B,^p - D^D^^p) = 
= e^°'^''n^{BaB„-i3 — DaD^-p 



(80) 



where the semi-colon stands fo r covariant differ entiation. 
The corresponding expression in iMcKinnevI (|2006l ) is a little 
bit different because it includes the term B" (np-^a+na-^p) , 
which equals to zero. Collecting all the results, we obtain 



.7" = ^e''^f''n^DpBs + 



B" 



B2 



(81) 
(82) 



It is easy to verify that in the 3-1-1 notation Ea l81l is identical 
to Eg 1701 which does not include neither the lapse function 
nor the shift vector, nor the time derivatives of B and D. 



4 ACKNOWLEDGMENTS 

It is a pleasure to thank Maxim Lyutikov and Jonathan 
McKinney for stimulating discussions. 



© 0000 RAS, MNRAS 000, 000-000 



